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The rolling motion of a rigid cylinder on an inclined flat viscous surface is investigated and the nonlinear 
resistance force against rolling, F R (y), is derived. For small velocities F R {v) increases with velocity due 
to increasing deformation rate of the surface material. For larger velocity it decreases with velocity due to 
decreasing contact area between the rolling cylinder and the deformed surface. The cylinder is, moreover, 
subjected to a viscous drag force and stochastic fluctuations due to a surrounding medium (air). For this 
system, in a wide range of parameters we observe bistability of the rolling motion. Depending on the 
material parameters, increasing the noise level may lead to increasing or decreasing average velocity. 



1 INTRODUCTION 

Rolling friction belongs to the most impor- 
tant dissipation mechanisms which may sensi- 
tively affect the dynamical behavior of granu- 
lar systems. Since 1785 when Vince described 
systematic experiments to determine the na- 
ture of friction laws , rolling friction has 
been investigated by many scientists due to 
its great importance in engineering and natu- 
ral sciences, e.g. (ICzichos 19781 IBarwell 19791 
LSzeri 19801 IDowson et al.T994l |Shpenkov 1995] 
Rab inowicz 1965I) . 

Experimentally it was observed that the rolling 
friction coefficient depends non-monotonically 
on the velocity (IQlson et al. 20031 TYi et al. 20021 
IGustafsson 19971 |Ray 1997[ IBakker et al. 19871 
IBurckhardt 1993l|Yi & Jeong 1998| ): For small ve- 
locities the rolling friction force increases with 
velocity, while for fast motion it decays. In 
the cas e of a soft cylinder ro lling on a hard 
plane (Brilli antov & Poschel 199 8). the contact 
surface between the bodies is flat which allows 



for the application of Hertz' contact theory ex- 
tended to the contact of viscoelastic particles 
(Brilliant ov et al. 1996b . In the opposite case, this 
assumption is not justified since the plane shape 
follows the shape of the rolling body in the contact 
area. Here we investigate rolling of a hard cylinder 
on a soft plane and bistability of rolling motion as 
a consequence of the resulting friction law. 

2 MODEL 

Consider a cylinder of radius R, mass M and 
moment of inertia I which rolls at velocity v on 
an inclined plane (quantities such as M, I, as 
well as the forces F R , F ex and similar quantities, 
are given as per unit length of the cylinder, 
i.e. as line densities), see Fig. d For certain 
materials it has been shown that surface effects 
such as adhesion may have significant influence 
on rolling friction (e.g. ( |Barquins et al. 19781 

I Deryaguin & Toporov 1994[ 

Kendall 19751 Fuller & Roberts 19811 

IRoberts & Tho mas 1975I0 . On the other hand, 



the contact area (x_,a; + ): 
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y(x) = R-h-y/B*-(x-x c ) 2 *i ( X 2 * e) -h. 

(2) 

where x c (t) is the position of the cylinder, h = 
—Vmin = —y (x c ) is the penetration depth of the 
cylinder and R^> \x — x c \ was assumed. Thus, in 
the stationary state, x c = vt, we obtain 



Figure 1 . A rigid cylinder rolls on an inclined plane built of 
uncoupled damped springs. The cylinder contacts the surface 
of the deformed plane in the region (£_,£+) which moves 
along with the cylinder. 

for viscoelastic materials it was reported that 
rolling friction is due very little to surface 
interactions, i.e., the major part is due to defor- 
mation losses within the bulk of the material 
(ITabor 19551 ITabor 19521 lAtack & Tabor 19581 
Drutowski 19 59). We will, therefore, neglect 
surface effects. The surface is modeled by non- 
interacting springs (f or the justification of this 
approximation see (Posc heTet al. 19991) ). The 
properties of the plane are described by the 
coefficients k, 7, and m: kd£ and 7d£ are the 
elastic and dissipative forces and md£ is the mass 
of the springs in the region d£. For X- < x < x + 
the surface is deformed by the cylinder. Assuming 
Hooke's law, the dynamics of the surface is 
described by 



my{x) + jy (x) + ky(x) = f (x,t) , (1) 



where / (x, t) is the density of the force which acts 
onto the springs in the region of contact. This force 
has to be provided by the cylinder. In a wide range 
of material parameters, the condition k < j 2 /4:k 2 
assures that the contact area is continuous, i.e., 
the cylinder contacts the surface everywhere in the 
interval ( (IPoschel et al. 199 9 ) for de- 

tails. 

From geometry we find the deformation y (x) in 



y(X) = -v- 



R 



y(x) 



and, therefore, 



/ U) = _^ _ lit + 
J 2iT R^ R 



R (3) 



hk . (4) 



with £ = x — x c . The boundary of the contact area 
at the front side of the cylinder in the direction of 
motion follows from geometry (in the co-moving 
frame), 

£+ = V2Rh (5) 

while the other boundary follows from the contact 
condition / (£_) = 0: 
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(6) 



The springs at x = x + are special: at time t — 5 
(5 — > +0) their velocity is zero, while an infinites- 
imal time later it is finite. This singularity may be 
attributed to a finite force F' N acting at the point 
x + . During the time dt the cylinder moves by vdt 
and accelerates springs of mass mvdt, therefore 



dp 
"dt 



— = — y [x + ) mv = \ —mv 




(7) 



The dissipation of energy due to the plane defor- 
mation and the instantaneous acceleration of the 
plane material at £ + in the instant of first contact 
reads 
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<%f{Z)y(t)--?v 2 (Z+)v = -vFj 



R ■ 



(8) 
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defining the force F R which acts against rolling 
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(9) 



To evaluate F R we have to determine the penetra- 
tion depth h(v). It is given implicitly by equilibrat- 
ing the total force exerted by the plane to the cylin- 
der and the weight of the cylinder, 



J f(0^ + F' N = Mgcosa. 



(10) 



where /(£) and F' N are given by Eqs. © and 
©. For an analytic evaluation of F R we as- 
sume = 0. This approximation means that 
the surface recovers so slowly, that the recover- 
ing surface does not transmit energy back to the 
cylinder! Gre enwood et al. 19611) . Since for realis- 
tic parameters (see below) the effects of interest, 
namely bistability of the rolling velocity and noise 
controlled velocity, occur at rather large velocity, 
this approximation is justified. Then Eq. (flOb re- 
duces to Mucosa = jvh + ^VMkh 3 / 2 , which, 
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for large damping, vy ^> kv2Rh, yields 



h(v) 



Mg cos a 



2RM g cos a 

-yv 



With the same arguments from Eq. 
the rolling friction force 



we obtain 



kh 2 



R 



mr2h + ^V2Rh. 
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(12) 



Figure |2] shows F R (v) as given by Eq. dT2b 
together with the numerical solution of Eqs. 
(I5I6I9I10I) for realistic material parameters, cor- 
responding to soft rubber (m =100kg/m 2 , k = 



10 7 kg/m 2 /s 2 
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5 x 10 & kg/m7s). The faster 



the cylinder moves, the more efficiently en- 
ergy is dissipated due to increasing deforma- 
tion rate. On the other hand, the faster it 
moves the smaller is the contact area, i.e. 
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Figure 2. Rolling friction force Fr(v). Full line: numeri- 
cal solution, dashed line: approximative theory, Eqs. il 11 and 
lO- 
the amount of deformed material decreases. 
Consequently, F R depends non-monotonically 
on the velocit y, which agrees with experi- 
mental results (IQlson et al. 20031 lYi et al. 20021 
IGustafsson 19971 |Ray 19971 IBakker et al. 19871 
IBurckhardt 19931 |yT& JeongJ998b. 



3 CYLINDER ROLLING DOWN AN IN- 
CLINE 

The rolling friction force F R is now applied to 
describe the motion of a cylinder rolling down a 
plane inclined by the angle a. The cylinder is sub- 
jected to an external driving force F ex = Mg sin a 
(g is the acceleration of gravity). The rolling fric- 
tion force F R and the viscous drag force F D due 
to the surrounding air counteract this motion. We 
also assume that the tangential force acting be- 
tween the cylinder and the surface at the contact 
area is strong enough to keep it from sliding. New- 
ton's equation for the cylinder reads 



(M + I/R) v = -F D (v) - F R (v) + F & 

The viscous drag force is given by 

F D = Av + Bv 2 . 



+ C(t). 

(13) 
(14) 



Detailed analysis dPoschel et al. 20051) suggests 
A = and B = 0.2kg/m 2 . Figure|3]shows the total 
force acting on the cylinder, F(v) = —Fd — F R + 
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Figure 3. Total force F(v) for different angle of inclina- 
tion, 1.5°, 2°, and 2.5°. Full lines: numerical solution, dashed 
lines: approximative theory. 

The steady state condition, v — 0, may be ful- 
filled either for only one velocity (top and bot- 
tom curves in Fig. |3j or for three different veloc- 
ities (middle curves). The former case implies a 
unique stationary velocity, while the latter allows 
for three stationary velocities. Only two of them, 
the smallest and the largest, correspond to stable 
motion (dF/dv < 0). Modifying the inclination a, 
the system may transit from one stable solution to 
the other. Figure|4]shows the bifurcation diagram. 

As shown in Figs. |2|4] the analytical theory 
(dashed curves) agrees well with the numerical re- 
sults if the damping parameter 7 is large. 




a 



Figure 4. Bifurcation diagram for stable velocity states (full 
line: numerical solution, dashed line: analytical approxima- 
tion). 



4 NOISE-INDUCED JUMPS 
Up to now, we did not consider noise, which is 
always present in realistic systems. The stochas- 
tic force ((t) describing fluctuations in the media 
is modeled by Gaussian white noise of zero av- 
erage, (((*)) = 0, and intensity a: (C(*)C(*0) = 
a 2 S(t — t'). If the system has only one stable ve- 
locity, the addition of noise does not change the 
motion of the cylinder qualitatively. In this case the 
velocity fluctuates around the average value given 
by the steady state condition v — 0. Figure |5| (top 
row) shows the velocity of the cylinder as obtained 
from numerical integration of the stochastic equa- 
tion (fT3t . The corresponding velocity distribution 
reveals a single peak (Fig. bottom row). 
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Figure 5. Top: velocity of the cylinder subjected to noise 
(<t =50kg 2 /s s ). Bottom: corresponding velocity distribu- 
tions. From left to right: a = 4.87°, 4.9°, 4.93°. The param- 
eters are 7 = 0.8 x 10 5 kg/m 2 /s, k = 10 6 kg/m 2 /s, B = 
0.65kg/m 2 . 

For parameters corresponding to bistable veloc- 
ity, the presence of noise changes the system quali- 
tatively due to stochastic jumps between the meta- 
stable velocities, Fig. |5] (middle column). Conse- 
quently, the velocity distribution has two well sep- 
arated peaks. 

Due to the nonlinear dependence of rolling fric- 
tion on the velocity, the average velocity may in- 
crease or decrease with increasing noise level, Fig. 
|6j This phenomenon may be understood by ana- 
lyzing the potential U(v) = —dF/dv. Depending 
on the parameters, U(v) has a double- or single- 
well shape dPoschel et al. 2005b . corresponding to 
the stationary state(s) of the velocity, see Fig.|3] In 
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Figure 6. Mean velocity of the cylinder over the noise inten- 
sity, (a) a = 4.87°, (b) a = 4.93°. 

both cases, U (v) is an asymmetric function in the 
surrounding of its minima where F(v) = 0. In the 
absence of noise the velocity remains in its (meta- 
)stable state. Subjected to noise, however, the ve- 
locity of the cylinder fluctuates around this mini- 
mum. If U (v) close to the minimum is steeper in 
the direction of lower velocities than in the direc- 
tion of higher velocities, the average velocity will 
be shifted towards higher velocities and, thus, in- 
creases with increasing noise level (Fig. |6k). In 
the opposite case, large noise impedes rolling (Fig. 

5 SUMMARY 

We studied the rolling motion of a hard cylinder 
on an inclined viscoelastic plane in the presence 
of a surrounding medium. For certain realistic pa- 
rameters the stationary velocity of the cylinder is 
bistable. For large damping of the plane's viscous 
deformation the numerical results agree well with 
theory. 

In the presence of noise as is unavoidable in 
any realistic experiment, by means of numerical 
simulations we found noise-induced transitions be- 
tween the meta- stable velocities. Depending on the 
system parameters, increasing noise level may ac- 
celerate or decelerate the rolling motion. 

The described effects may be important for tech- 
nical systems where the presence of noise may lead 
to an effective increase of the mobility of a rolling 
body which is driven by an external force. 
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